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Abstract 

In this work, we calculate energy loss of heavy quark (charm and bottom) due 
to elastic collisions and gluon radiation in hot/dense medium. The collisional 
energy loss has been obtained using QCD calculations. The radiative energy 
loss is calculated using reaction operator formalism and generalized dead cone 
approach. We rederive the energy loss expression using same assumptions as 
generalized dead cone approach but obtain slightly different results. We also 
improve the model employed to calculate path length and the system evolution. 
The nuclear modification factors Raa including shadowing and energy loss are 
evaluated for B and D mesons and are compared with the measurements in 
PbPb collision at yjS nn = 2.76 TeV and with the D meson and Heavy flavour 
(HF) electrons measurements in AuAu collision at yjS nn = 200 GeV. The radia¬ 
tive energy loss calculated by reaction operator formalism added with collisional 
energy loss describes the RHIC HF electron suppression in high pt range. It 
also describes the LHC measurement of B meson suppression but overestimates 
the suppression of D meson. The radiative energy loss from generalized dead 
cone approach describes the charm suppression at both RHIC as well as LHC 
energies and requires energy loss due to collisions to be added in order to de¬ 
scribe the bottom suppression at LHC. 
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1. Introduction 

The heavy ion collisions at ultra relativistic energy create matter with high 
energy density required to form Quark Gluon Plasma (QGP). Relativistic Heavy 
Ion Collider (RHIC) and Large Hadron Collider (LHC) are designed to create 
and explore QGP. Many measurements at RHIC and LHC already point to the 
formation of QGP pi. The heavy quarks (charm and bottom) are produced in 
hard partonic interactions in heavy ion collisions and their initial momentum 
distribution can be calculated from pQCD [2]. While traversing the hot/dense 
medium formed in the collisions, these quarks loose energy either due to the 
elastic collisions with the plasma constituents or by radiating a gluon or both. 
There are several formulations to calculate collisional 0 n 0 0 m as well as 
radiative energy loss 0 ® eh inj. For a review of many of these formalism 
see Ref. [Hj. At high parton energies, the radiative energy loss becomes much 
larger than the collisional energy loss but at lower energies, these two processes 
can contribute equally with the collisional energy loss being the dominant for 
small values of the parton energy m 

There are many heavy quark energy loss models, each having specific set of 
simplifications/assumptions. The model by Gyulassy, Levai and Vitev (GLV) 
pn m is based on a systematic expansion of the energy loss in terms of the 
number of scatterings and generally leading order term is included in the cur¬ 
rent calculations. The medium is characterized by two parameters, the density 
of scattering centers or mean free path and Debye screening mass. Such an 
approach includes the interference between vacuum and medium induced radi¬ 
ation. This formalism was then extended to obtain the energy loss for heavy 
quarks in Ref. m and was simplified for the first order of opacity expansion in 
Ref. pi] 


2 



In this work, we calculate the radiative energy loss of heavy quarks (both 
charm or bottom quark) using reaction operator formalism DGLV (Djordjevic, 
Gyulassy, Levai and Vitev) [TUI HT1115] and using generalized dead cone approach 
AJMS (Abir, Jamil, Mustafa and Srivastava) [T7J. We rederive the energy loss 
expression using same assumptions as generalized dead cone approach but obtain 
slightly different results. We also improve the model employed to calculate path 
length and the system evolution. The collisional energy loss has been calculated 
using Peigne and Peshier formalism [Jj. The nuclear modification factors Raa 
including shadowing and energy loss are evaluated for B and D mesons and are 
compared with the measurements in PbPb collision at a/s]v]v = 2.76 TeV and 
with the HF electron measurement of PHENIX and D meson measurements of 
STAR in AuAu collision at yj snn = 200 GeV. 

2. Heavy Quark Production by Hard Processes 

The production cross sections of cc and bb pairs are calculated to NLO in 
pQCD using the CT10 parton densities [T5] . We use the same set of parameters 
as that of Ref. m which are obtained by fitting the energy dependence of 
open heavy flavor production to the measured total cross sections. The charm 
quark mass and scale parameters used are m c = 1.27 GeV, Hp/mpc = 2.10, and 
Hr/hitc = 1-60 [T9] , The bottom quark mass and scale parameters are mb = 
4.65 GeV, pp/mpb = 1-40, and / iR/mpb = 1.10. Here hf is the factorization 
scale, hr is the renormalization scale and mp = \JM 2 + pi r . The central EPS09 
NLO parameter set |20i is used to calculate the modifications of the parton 
distribution functions (nPDF) in heavy ion collisions, referred as shadowing 
effects. 

For the fragmentation of heavy quarks into mesons, Peterson fragmentation 
function is used which is given as follows [HQ 

N 

Do(2) = *[1 - (IM - « a /d - *)r (1) 

Here z = p® /Pp and N is normalization constant which is fixed by summing 
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over all hadrons containing heavy quark, 

J dzD Q ( z ) = !• ( 2 ) 

We take e c = 0.016 and = 0.0012. The schemes of D meson to electron decay 
(BR = 10.3 %) and B meson to J/^> decay (BR = 1.1 %) are obtained by Pythia 
simulations f?21 . 

Figure [T] shows the NLO calculations of differential cross section of single 
electrons from D mesons as a function of the transverse momentum px in pp 
collision at y/s = 200 GeV compared with the PHENIX measurements of single 
electrons from heavy flavour [23]. As shown in the Data/Theory panel the 
agreement between the data and the calculations is not very good but since at 
Pt above 2 GeV/c the shapes of the calculations and the data are same, it does 
not affect the Raa calculations due to energy loss. 

Figure [ 2 ] shows the NLO calculations of differential cross section of D° 
mesons as a function of the transverse momentum px in pp collision at y/s = 2.76 
TeV compared with ALICE measurements of D° and D + mesons [23]. Here also 
Data/Theory panel shows the agreement between the data and the calculations 
is not very good but since at px above 2 GeV/c the shapes of the calculations 
and the data are same, it does not affect the Raa calculations due to energy 
loss. 

Figure [3] shows the NLO calculations of differential cross section of inclusive 
J/i/j coming from B mesons as a function of the transverse momentum px in pp 
collision at y/s = 2.76 TeV. 

3. Collisional Energy Loss 

The QCD calculation of the rate of energy loss of heavy quark per unit 
distance ( dE/dx ) in QGP is given by Braaten and Thoma [3]. Their formalism 
is an extension of QED calculation of dE/dx for a muon |T| which assumes that 
the momentum exchange q <C E. Such an assumption is not valid in the domain 
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Figure 1: (color online): The NLO calculations of differential cross section of single electrons 
from D mesons as a function of the transverse momentum pt in pp collision at >/s = 200 
GeV. The data is from PHENIX measurements of single electrons from heavy flavour [23]. 



p r (GeV/c) 


Figure 2: (color online): The NLO calculations of differential cross section of D° mesons as a 
function of the transverse momentum pt in pp collision at y/s = 2.76 TeV. The data is from 
ALICE measurements of D° and mesons m- 
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Figure 3: (color online): The NLO calculations of differential cross section of inclusive J/ip 
coming from B mesons as a function of the transverse momentum pt in pp collision at y/s = 
2.76 TeV. 


when the energy of the heavy quark E S> M 2 /T, where M is the mass of the 
heavy quark. Peigne and Peshier (PP) !7j extended this calculation which is 
valid in the domain E M 2 /T to give the expression for dE/dx as 


dE 

dx 


4:Tra 2 T 2 




log 


^ V a ' 


2, ET 
9 ° g M 2 


' c (Nf) 


(3) 


Here fi g = \nasT 2 (l + Nf/6^J is the Debye screening mass and c(Nf) ~ 
0.1467V/ + 0.05 . a s (= 0.3) is the fine structure splitting constant for strong 
interaction and TV/ is the number of quark flavours. 


4. Radiative Energy Loss 
4-1. DGLV Formalism 

The energy loss of fast partons is dominated by radiation of gluons. The reaction 
operator formalism is used in Ref. m to obtain the energy loss due to gluon 
radiation for light quark jets. Analytical expression is obtained for energy loss in 
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powers of gluon opacity (L/ A) where A is the mean free path of the quark and L 
is the path length traversed in the medium. This formalism was then extended 
to obtain the energy loss for heavy quarks in Ref. m and was simplified for 
the first order of opacity expansion in Ref. m- The expression of the average 
radiative energy loss of heavy quark is given in appendix A. 


4-2. Generalized Dead Cone Approach (AJMS) 

The rate of radiative energy loss of a heavy quark with energy E due to the 
inelastic scattering with the medium is calculated as 


dE < oj > 
dx A 

where < w > is the mean energy of the emitted gluons. 

The probability of gluon emission off a heavy quark is written as HZI 


(4) 


drjg ‘ZCacxs T) 

dr) du 7 r ui 

where Ca{= 3) is the Casimir operator in QCD and w is related to the transverse 
momentum of the emitted gluons k± by the relation k± = to sin 9, where 6 is 
the emission angle. T> is the generalised dead cone which can be written as [251 


V = ^1 + e 2v ^j , r)=— In tan . (6) 

Here s is mandelstam variable which is related to the energy E and mass M of 
heavy quark by the relation, s = 2 E 2 + 2 E\/E 2 — M 2 — M 2 . 

The mean energy of the emitted gluon can be written as m 

^_fSt ujdr ld^ JdufVdr) 

" JSbdndu, Ib^fVdr)- U 

The mean free path length A is calculated as m 


Pq a Qq(q)^>-Qq(q)g + Pg a Qg^Qgg i 

(8) 

, 9 . 

\Pq + ^ Pg) a 2^3, 

(9) 

PQGP &2^>3 ■ 

(10) 
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The total cross section of the process 2 —► 3 is calculated as [22] 


(T2->3 — 4 Ca a? 


1 


2J2 


1 


duj 


V dr] . 


(ID 


_ {dl ) 2 ~~ J w 

Here q± is the transverse momentum of the exchanged gluon. Using Eqs. 
0- §. 0 and ( fTlj ) and assigning the limits of the variables of q 2 , ui and 
r] we get 

dE nt 3 [ q ±\•»“ 1 

— = 24 a s pQ GP / 






2 ? dr] . ( 12 ) 


- / / 2 \2 

“ I min \'i±) J Tjmin 

Here we have put Ca= 3 and factor of 2 is used to cover both upper and lower 
hemispheres of 77 . 


Equation (12 1 is solved to get following result (see details in appendix B) which 
we call corrected AJMS result 


dE , 1 , , 

— = 24 a s PQGP-[1- 




Here 


T(6) = 25-~ log 


1 + AU e 25 
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sinh(2<5) 


M 2 e _2A 


1 + 2 
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<5 = 2 lo S 


(1-/30 l0S (/3i 


I + 4 /I- 


(l~/3i) 
lQ g(^) , 


(13) 


(14) 


(15) 


The above results differs with the original AJMS calculation m where the 
T(5) term is given by 


m = 


5 = 


1 / 1 + Ml e 25 ' 


Yr cosh ((5) 


l + 2f cosh (J) + ^ 
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5. Model For QGP Evolution 


To estimate the energy loss suffered by the heavy quark, it is crucial to 
calculate its path length which it travel in the medium. Let us assume that 
the heavy quark is produced at a point (r, <f>) in heavy ion collision, moves at 
an angle <f> with respect to r in the transverse plane. If R is the radius of the 
colliding nuclei, then the distance d covered by the heavy quark in the plasma 
is given m by 



(17) 


The average distance travelled by the heavy quark in the plasma 


L = J 0 * JcT d(^),r) p{\r\) p(\r- b\) r dr d</> 

f 0 R /(T p( i^i) ~ &l) r dr ^ 


(18) 


Here p(|f|) is the density of nucleus assumed to be a sharp sphere with radius 


R = 1.1 d 1 / 3 . If the velocity of the heavy quark is vt = Pt/wit, where tut is 
the transverse mass, the effective path length L e // is obtained as 


(19) 


L e ff = min L , vt x 17 


The evolution of the system for each centrality bin is governed by an isen- 
tropic cylindrical expansion with prescription given in Ref. The entropy 

conservation condition s(T) V(t) = s(Tq) V(ro) and equation of state obtained 
by Lattice QCD along with hadronic resonance are used to obtain temperature 
as a function of proper time [29] . The transverse size R for a given centrality with 
number of participant N part is obtained as R(N part ) = Ra \J 2 A/N part , where 
Ra is radius of the nucleus. The initial entropy density s(r 0 ) is 



( 20 ) 


Here a m = 5 is a constant which relates the total entropy with the multiplic¬ 


ity [30). The initial volume V(t 0 ) = n [i?(lV part )] 2 r 0 and measured values of 


dN/dr 1 for LHC J2I) and for RHIC [32! are use d for a given centrality. 
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Table 1: Parameters of QGP evolution model 


Model 

Present 

Present 

Present 

Model 

mm3 

Model 

mm3 

Centrality (%) 

0-10 

0-100 

0-20 

0-10 

0-20 

Experiment 

RHIC 

LHC 

LHC 

RHIC 

LHC 

b (fm) 

3.26 

9.72 

4.70 

0.0 

0.0 

-^part 

329 

113 

308 

- 

- 

dN/dr] 

623 

360 

1206 

- 

- 

L( fm) 

5.63 

4.3 

5.62 

5.78 

6.14 

T 0 (fm/c) 

0.6 

0.3 

0.3 

0.2 

0.2 

T f (fm/c) 

3.0 

6.0 

6.0 

2.63 

5.90 

T 0 (GeV) 

0.303 

0.450 

0.481 

0.400 

0.525 


We calculate the energy loss which is then averaged over the temperature 
evolution. Various parameters used and calculated in our model for different 
centralities such as number of participant iVp art , measured dN/drj, calculated 
average path length (T), initial time (to), QGP life time (tj) and initial tem¬ 
perature (To) are given in Table [lj The parameters used in the earlier model 
mum are also given. 

6. Results and Discussions 

Figure [4] shows energy loss of charm quark as a function of energy of quark 
for AuAu collision at y/sNN =200 GeV using PP, DGLV, AJMS and corrected 
AJMS formalisms. Figure [5] is the same for bottom quark. It can be seen that 
the collisional energy loss is similar in magnitude for charm and bottom quark. 
The radiative energy loss calculated by AJMS is larger than that by DGLV. 
This difference is more pronounced for bottom quark. Radiative energy loss 
of bottom quark by AJMS starts dominating collisional energy loss at quark 
energy above 11 GeV whereas the DGLV energy loss remains below collisional 
energy loss upto 25 GeV of bottom quark energy. Figure [6] shows energy loss of 
charm quark as a function of energy of quark for PbPb collision at y/sNN= 2-76 
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TeV using PP, DGLV, AJMS and corrected AJMS formalisms. Figure [7] is the 
same for bottom quark. When we move from RHIC to LHC both the collisional 
as well as radiative energy loss increase. The radiative energy loss of charm 
quark calculated by AJMS and DGLV are similar in magnitudes but the DGLV 
energy loss increases more steeply with quark energy. For bottom quark, the 
AJMS energy loss is much larger than the DGLV energy loss. Radiative energy 
loss of bottom quark by AJMS starts dominating collisional energy loss at quark 
energy above 10 GeV whereas the DGLV energy loss remains below collisional 
energy loss upto 22 GeV of bottom quark energy. 



Figure 4: (color online): The energy loss dE/dx as a function of energy of charm quark ob¬ 
tained using PP, DGLV, AJMS and corrected AJMS (Present) calculations for AuAu collision 
at y/sjjN=200 GeV. 


Figure [8] shows nuclear modification factor Raa of single electrons from D 
meson as a function of transverse momentum px obtained using energy loss 
(DGLV, AJMS, corrected AJMS and DGLV+PP calculations) and shadowing 
in AuAu collision at v /s)v)v=200 GeV. The data is PHENIX measurements 
of heavy flavour (HF) electrons S3] . We observe that radiative energy loss by 
DGLV added to collisional energy loss by PP describes the PHENIX data at high 
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Figure 5: (color online): The energy loss cLE/dx as a function of energy of bottom quark ob¬ 
tained using PP, DGLV, AJMS and corrected AJMS (Present) calculations for AuAu collision 
at y/sjjN=200 GeV. 



Figure 6: (color online): The energy loss dE/dx as a function of energy of charm quark ob¬ 
tained using PP, DGLV, AJMS and corrected AJMS (Present) calculations for PbPb collision 
at y/sNN= 2.76 TeV. 


12 







Figure 7: (color online): The energy loss dE/dx as a function of energy of bottom quark 
obtained using PP, DGLV, AJMS and corrected AJMS (Present) calculations in PbPb collision 
at - v /sjvat=2.76 TeV. 


Pt range. The radiative energy loss by AJMS and corrected AJMS reproduce 
the data without adding energy loss due to collisions. 

Figure [9] shows the nuclear modification factor Raa of D meson as a function 
of transverse momentum pt obtained using energy loss energy loss (DGLV, 
AJMS, corrected AJMS and DGLV+PP calculations) and shadowing in AuAu 
collision at ■ s /snn =200 GeV. The data is STAR measurements of D° mesons 
[33;. We observe that the radiative energy loss by AJMS and corrected AJMS 
reproduce the data without adding energy loss due to collisions. The radiative 
energy loss by DGLV added to collisional energy loss by PP describes the STAR 
data at high p t range. 


Figure 10 shows the nuclear modification factor Raa of D° mesons as a 
function of transverse momentum obtained using radiative energy loss (cor¬ 
rected AJMS calculations) calculated with old and new evolution models and 
shadowing in PbPb collision at , /snn= 2.76 TeV. The data is ALICE measure¬ 
ments of D° mesons [31] The value of Raa depends on energy loss model as 
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well as evolution model used to calculate the pathlength. 

shows the nuclear modification factor Raa of D Q mesons as a 


Figure 
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function of transverse momentum px obtained using energy loss (DGLV, AJMS, 
corrected AJMS and PP+DGLV calculations) and shadowing in PbPb collision 
at ■ s /snn= 2.76 TeV. The data is ALICE measurements of D° mesons i.Tl . 
AJMS , corrected AJMS and DGLV calculations produce similar suppression in 
high px range. When we add radiative and collisional energy loss (PP+DGLV) 
it overestimates the measured suppression of D meson. 


Figure 12 shows the nuclear modification factor Raa inclusive J/i/j coming 
from B mesons as a function of transverse momentum px obtained using energy 
loss (DGLV, AJMS, corrected AJMS and PP+DGLV calculations) and shad¬ 
owing in PbPb collision at v /satjv= 2.76 TeV. The data is CMS measurements 
of J/ip mesons from B decays (35;. We observe that sum of radiative energy 
loss (DGLV) and collisions energy loss (PP) underestimates the B meson sup¬ 
pression. The sum of radiative energy loss by corrected AJMS and collisions 
energy loss slightly overestimates the suppression. More accurate data in larger 
px range would help distinguish the models more clearly. 

Figure [13] shows the nuclear modification factor Raa of single electrons from 
D meson as a function of the number of participant N part obtained using en¬ 
ergy loss (DGLV, AJMS, corrected AJMS (Present) and PP+DGLV calcula¬ 
tions) and shadowing in AuAu collision at y/sNN =200 GeV compared with the 
PHENIX measurements of heavy flavour (HF) electrons [23]. The radiative en¬ 
ergy loss by DGLV added to collisional energy loss by PP slightly overestimates 
the suppression. The radiative energy loss by corrected AJMS describes the 
data without adding energy loss due to collisions. 

Figure [14] shows the nuclear modification factor Raa of single electrons from 
D meson as a function of the number of participant N part obtained using en¬ 
ergy loss (DGLV, AJMS, corrected AJMS and PP+DGLV calculations) and 
shadowing in AuAu collision at snn = 200 GeV compared with the STAR 
measurements of D mesons [33]. We observe that the radiative energy loss by 
AJMS and corrected AJMS describe the data without adding energy loss due 
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Figure 8: (color online): Nuclear modification factor Raa of single electrons from D meson as 
a function of transverse momentum pt obtained using energy loss (DGLV, AJMS, corrected 
AJMS (Present) and DGLV+PP calculations) and shadowing in AuAu collision at ^/SjViV—200 
GeV. The data is from PHENIX measurements of heavy flavour (HF) electrons m- 
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P T (GeV) 


Figure 9: (color online): Nuclear modification factor Raa of D meson as a function of trans¬ 
verse momentum pjp obtained using energy loss (DGLV, AJMS, corrected AJMS (Present) 
and DGLV+PP calculations) and shadowing in AuAu collision at y / sW=200 GeV. The data 
is from STAR measurements of D mesons m- 


15 






















2 » 1 1 1 I 1 ^3 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 

- PbPb = 2.76 TeV 0-20% centrality |y| < 0.5 


1.8 

1 . 6 ; 

1.4; 

1 . 2 ; 

< I 

< 1 
c ■- 

0 . 8 : 




0.4 


0.2 


• D° meson (ALICE) 
| | Syst. Unc. 

Shadowing Only 
Shadowing + Eloss + OldEvol. 

Shadowing + Eloss + Present_Evol. 




8 10 12 14 16 18 20 

P T (GeV) 


Figure 10: (color online): Nuclear modification factor Raa °f mesons as a function 
of transverse momentum obtained using radiative energy loss (corrected AJMS (Present) 
calculations) calculated with old and new evolution models and shadowing in PbPb collision 
at y/sNN= 2.76 TeV. The data is from ALICE measurements of D° mesons m- 


2p 
1.8 — 
1.6 
1.4; 

1.2 T 

< 

< I — 

0C 1 
0 . 8 - 

0.6 
0.4 
0.2 r 

°2 


PbPb \ls^ = 2.76 TeV 0-20% centrality |y| < 0.5 

• D° meson (ALICES 

- Shadowing Only ,-. 

. Shadowing + DGLV I-1 Syst. Unc. 

- Shadowing + AJMS 

- Shadowing + Present 

Shadowing + PP + DGLV 



111 ■ 11111111 ■ 111111111 ■ 111111 n 

6 8 10 12 14 16 18 20 

P T (GeV) 


Figure 11: (color online): Nuclear modification factor Raa °f mesons as a function of trans¬ 
verse momentum pt obtained using energy loss (DGLV, AJMS, corrected AJMS (Present) 
and PP+DGLV calculations) and shadowing in PbPb collision at y/SNN= 2.76 TeV. The data 
is from ALICE measurements of D° mesons m- 
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Figure 12: (color online): Nuclear modification factor Raa inclusive J/'ip coming from B 
mesons as a function of transverse momentum pt obtained using energy loss (DGLV, AJMS, 
corrected AJMS (Present) and PP+DGLV calculations) and shadowing in PbPb collision at 
yjs nn = 2.76 TeV. The data is from CMS measurements of J/'ip mesons from B decays ED. 


to collisions. The energy loss by DGLV does not describe the data. 


Figure 15 shows the the nuclear modification factor Raa of D° mesons as a 
function of the number of participant N part obtained using energy loss (DGLV, 
AJMS, corrected AJMS (Present) and PP+DGLV calculations) and shadowing 
in the PbPb collision at i/satjv=2.76 TeV compared with ALICE measurements 
of D° mesons [34[. We observe that the radiative energy loss by DGLV, AJMS 
and corrected AJMS describe the ALICE data. The radiative energy loss by 
DGLV added to collisional energy loss by PP overestimates the D° suppression. 


Figure 16 shows the nuclear modification factor Raa of inclusive J/ij) com¬ 
ing from B mesons as a function of the number of participant N part obtained 
using energy loss (DGLV, AJMS, corrected AJMS (Present), PP+DGLV and 
PP+corrected AJMS calculations) and shadowing in the PbPb collision at 
yjs nn = 2.76 TeV compared with the CMS measurements of J/i/’ mesons from 
B decays [35;. We observe that the radiative energy loss by DGLV added to 
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collisional energy loss by PP describes the CMS data very well. The radiative 
energy loss by corrected AJMS underestimates the suppression but with col¬ 
lisional energy loss added it overestimates the suppression. Both the models 
favour that B meson loose energy by radiation as well as collision processes. 



Figure 13: (color online): Nuclear modification factor Raa °f single electrons from D me¬ 
son as a function of the number of participant N par t obtained using energy loss (DGLV, 
AJMS, corrected AJMS (Present) and PP+DGLV calculations) and shadowing in AuAu col¬ 
lision at y/s nn=200 GeV. The data is from PHENIX measurements of heavy flavour (HF) 
electrons m- 


7. Conclusion 

We study the energy loss of heavy quark (charm and bottom) due to elastic 
collisions and gluon radiation in hot/dense medium. Results of Radiative energy 
loss obtained from two different formalisms namely DGLV and AJMS have been 
compared. The energy loss calculated by AJMS exceeds DGLV results. The col¬ 
lisional energy loss has been calculated using Peigne and Peshier formalism and 
is found to be similar in magnitude for charm and bottom quarks. The nuclear 
modification factors Raa including shadowing and energy loss are evaluated for 
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Figure 14: (color online): Nuclear modification factor Raa °f single electrons from D meson 
as a function of the number of participant N par t obtained using energy loss (DGLV, AJMS, 
corrected AJMS (Present) and PP+DGLV calculations) and shadowing in AuAu collision at 
y/s at tv— 200 GeV. The data is from STAR measurements of D mesons m- 



Figure 15: (color online): Nuclear modification factor Raa °f D° mesons as a function of 
the number of participant N par t obtained using energy loss (DGLV, AJMS, corrected AJMS 
(Present) and PP+DGLV calculations) and shadowing in the PbPb collision at ^/SiViV—2.76 
TeV. The data is from ALICE measurements of D° mesons m- 
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Figure 16: (color online): Nuclear modification factor Raa °f inclusive JxJj coming from B 
mesons as a function of the number of participant N par t obtained using energy loss (DGLV, 
AJMS, corrected AJMS (Present), PP+DGLV and PP+corrected AJMS calculations) and 
shadowing in the PbPb collision at y/s nn= 2.76 TeV. The data is from CMS measurements 
of J/'ip mesons from B decays ED- 
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B and D mesons and are compared with the measurements in PbPb collision at 
y/sNN = 2.76 TeV and with the HF electrons and D° meson measurements in 
AuAu collision at ^/sjq n = 200 GeV. The radiative plus collisional energy loss 
(PP+DGLV) describes the RHIC HF electron suppression in high px range. It 
also describes the LHC measurement of B meson suppression but overestimates 
the suppression of D meson. The radiative energy loss from generalized dead 
cone approach describes the charm suppression at both RHIC as well as LHC 
energies without requiring collisional energy loss. Both collision as well as ra¬ 
diative energy loss are required to explain the B meson suppression at LHC. 
Upcoming high luminosity PbPb collisions at ^/s^n = 5 TeV are expected to 
measure the heavy quarks in wider kinematic ranges which will provide much 
improved constrains for the processes of energy loss and models. 
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Appendix A: DGLV Formalism 
The average radiative energy loss of heavy quarks as 


A E 
~L~ 



where 


II 

CM 

m 2 g (l - x) + M 2 X 2 , A 1 = PgVQg + PqCTQq , 

(22) 

Pa = 

, 1.202 „ r , 1.202 

16 T 3 -2-, p q —■ 9 Nj T 3 -2-, 

7 T z 7 T z 

(23) 

°Qa = 

9ira 2 4 

2 a “' ' 

(24) 


Here Cf{= 4/3) determines the coupling strength of gluon to the massive quark 
with momentum p. p g and p q are the densities of gluons and quarks and m g = 
Pg/y/2 is the transverse gluon mass. 
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The function A, 

£? and 

C are given as follows 



A 

2/3 

ft 

V + <? 2 ), 


(25) 


B 


(/3 2 + K)(f3 2 Q~ + 

Q+Q+ + Q+ft) 

(0R\ 


/3 2 

(/3 2 (Q)1 - K) - Qj, 

A + Qg Q/t + fpfpj 



C 

2g- 

[/3 2 Ms(V ^ M 2 ) + /3 2 (/3 2 ^ fi)K 




+Q+(^ 4 - 2g 2 Q+) + / M (/3 2 (-/3 2 - 3g 2 + M 2 ) 




+2? 2 Q^ + 3/3 2 q"Q k 


(27) 

Here 








I< 

famax 2px(l 

- x), 

(28) 



Qi 

= q 2 ± Mg! <9^ 

— „ 2 + u 2 

- 1 - ^maxi 

(29) 



ft 

= fW,Q»,Q+), 

fn = /(Mgi Qk ’ Qk )> 

(30) 


/( 

x,y,z) 

= \/x A + 2x 2 y + 

Z 2 . 

(31) 


Appendix B: Corrected AJMS 


The integration of Eq. (12) are obtained as follows. 

The minimum values of q]_, lo and k 2 ^ are given by infra-red cut-off [T7I 


2 | ^ 2 ^ 7,2 I 

y_L I min ~ w min ~ K jJmin 

The maximum value of is calculated as 

<7i_ | max = C E T , 


fa 


(32) 


(33) 


where 


3 M 2 M 4 rM 2 + 6ET(l + ft)l 

2 ~~ 4 E T + 48 E 2 T 2 /3 0 ° S l-M 2 + 6 ET (1 -/3 0 )- 


and 



(34) 


( 35 ) 
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The maximum value of w is obtained as | E9j 


^max ^ 9_L ' 


( 36 ) 


The average of square of the transverse momentum q± is given in reference 

[Ml Eli as 

n 2 \ - n 2 \ rn 2 I 

(37) 


2 Q± | min 9_L I max 

< 9j_ >= 721-- lQ g 


911 max 


L 9ii 


9il max -9il min 

Putting q\\min from Eq. (32) and q±\ m ax from Eq. (33) in Eq. (37) 


2 <4 , 

< 91 >= 77 —log 


( 1 -A) L/ 3 X J ’ 


(38) 


where /3i = /i 2 /(C 7? T). Using the relation uj = k± cosh^, the Hnite cut on 
ui and k± leads to an inequality 


t^max 
k±\ min 


< cosh 77 . 


(39) 


The integration limits of 77 are calculated from Eq. (|32|), (|36|) and (|39|) as 

\v\ < log 


'<g± > , J <q1 > 


f4 




(40) 


We can write it as |r/| < 5, where S is obtained using equation (38) and (40) 

2 


5 =2 l0g 


1 


106 (£) \ 1 + f- 


(1 — Pi) 


(i-A) - va; cv* log(^) 

We can write the minimum and maximum value of 77 as 

Vmin = 'Hmax = 


(41) 


(42) 


Now we calculate the integrals in Eq. (12) which can be written as 


dE 

— = 24 a 3 s pqgp h h h ■ 


The first integration I\ is calculated as 

r G E T 


h = 


1 

A) 


(9 ^ 2 


dql = i (l - ft) ■ 

Pg 


(43) 


(44) 
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The second integration / 2 is calculated as 


h = 


did = \lg 


I 1 

(1 - Pi) 


o 


-1 . 


The third integration I 3 is calculated as 


I 3 = 


V dr] = 


1 


1-5 


(i + + 


dii 


1 + 


M 2 


log(e 4lS ) + log 


e 2 V 


log 


Ml e 2v 


1 + 


M 2 


e 2 V 


-5 


Ml a -25' 


1 +^~ e 


M_ ( e 25 _ e -25) 


1 + Ml e 25 y ^1 + Ml ( e 25 + e -2«) + ^ 

The integration / 3 is denoted as T{5) given by 


1 


F[S) = 25-- log 
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1 + Ml e 25 


M 


sinh(2<5) 


Ml e -25 


1 + 2 


M 2 


cosh(2d) 


M 4 


(45) 

(46) 

(47) 

^•) 

(49) 
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